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How to determine the symmetry of the superconducting order parameter is one of the important issues in novel 
superconductors, which include charge conjugated organic superconductors. We have proposed that the angular 
dependence of the thermal conductivity in a planar magnetic field provides a new window to look at the symmetry 
of the order parameter. After a brief summary of the quasiclassical approach we describe how the symmetry of 
the superconducting order parameter in Sr2Ru04, CeCoIns and k-(ET)2Cu(NCS)2 is determined. Also in some 
of experiments the phononic thermal conductivity plays the crucial role. 
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1. Introduction 

Since the discovery of the first organic super- 
conductor in (TMTSF) 2 PF 6 (or Bechgaard salts) 
PP, the symmetry in superconducting order pa- 
rameter has been one of the important issues |2I3| . 
Indeed the symmetry of the superconducting or- 
der parameter becomes one of the central issues 
after the establishment of d-wave symmetry in 
both hole and electron doped high T c supercon- 
ductors 4 5 6 7. Most likely the superconductiv- 
ity in (TMTSF) 2 X with X = C10 4 , PF 6 , etc. is 
of p-wave [S]. In particular, a flat Knight shift 
seen in a recent NMR experiment [§] is consis- 
tent with this picture. The remaining question is 
whether the p-wave superconductor belongs to ID 
representation [3] or 2D representation [Sj. After 
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a long controversy ^UJi d-w&ve superconductiv- 
ity is emerging in k-(ET) 2 salts |llll2ll3ll4IT5] . 
Here now the question is whether the symme- 
try is of d xy -wave as suggested by the theoret- 
ical works [16117118] or d^.^-wave as the re- 
cent STM study suggests |T2]- We shall give a 
somewhat surprising answer on this based on the 
recent angular dependent thermal conductivity 
data of k-(ET) 2 Cu(NCS) 2 by Izawa et al |T3]. 
The quasiparticle spectrum of all these new su- 
perconductors is well described by the BCS the- 
ory for nodal (or unconventional) superconduc- 
tors |6l7l2()j . In particular, there arc nodal excita- 
tions (i.e. the quasiparticles which inhabit in the 
vicinity of the nodal lines) which persist to low 
temperatures (i.e.T<< A where A is the super- 
conducting order parameter) . In the vortex state 
the quasiparticle spectrum is modified due to the 
supercurrent circling around individual vortices. 
In order to describe the quasiparticle spectrum 
in the vortex state, Volovik |2J has introduced a 
very simple method to evaluate the effect of the 
supercurrent within the quasiclassical approxima- 
tion. In particular, he has shown that the specific 
heat in the vortex state in nodal superconductors 
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is proportional to \fTl for H/H c2 « 1 where H 
is the magnetic field strength. This yH depen- 
dence has been seen in YBCO |22l23l24| . LSCO 
[25] . k-(ET) 2 salts Q2] and Sr 2 Ru0 4 [2^l27| . We 
shall show later that in the presence of impurity 
the above dependence may change as 
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for Sr2Ru04 in H// a-b. Similarly for d+s-wave 
in H// a-b, 
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Here iV(0) is the residual density of states in the 
presence of the impurity scattering and Nq is the 
one for the normal state. Also 7jv is the Som- 
mcrfcld coefficient, v — ^/v a v c and v a and v c are 
the Fermi velocities in the a-b plane and paral- 
lel to the c-axis, respectively. For d + s-wave 
we took A(k) oc cos(2</>) — r. 8 is the angle 
H makes from the a axis. Here x = v • q/A 
and v • q is called the Doppler shift, where v is 
the Fermi velocity and 2q is the pair momentum 
and < ... > means space average over the vor- 
tex lattice and v average over the nodal lines. 
Some of these details are given in |29|. In the 
above derivation, we have assumed that the sys- 
tem is in the clean limit ( e, T « VTA ) while 
Volovik's result applies for the superclean limit 
(\/fA < e, T). Here e = \v\feH is the char- 
acteristic magnetic energy. Also when we put r 
= in Eq.(2), we obtain the usual expression for 
d-w&ve superconductors. The specific heat in the 
clean limit in LSCO in H || c has been reported 
recently The quasiclassical approximation is 
extended to calculate the thermal conductivity 
in the vortex state |28I31I32I33I53| . Indeed we 
can now describe the angular dependent thermal 
conductivity observed in single crystals of YBCO 
138) consistently if we assume that the 



system is in the superclean limit and T>> e |34|. 
In the following we shall first review the theory 







Figure 1. |A(k)|'s for s-wave and some 
of unconventional superconductors are shown: 
(a) s-wave superconductor, A(k) = A (b) 
d a; 2_ !/ 2-state, A(k) ~ cos(2cf>) (c) 2D /-wave, 
A(k) - cos(cfc z )e ± ^ (d) £ 2 u-state, A(k) ~ 



cos t) sin 2 9e ±lcf 



limiting ourselves to the clean limit (e << %/AT). 
As to the result for the superclean limit read- 
ers may consult |27I34I39| . Also we consider the 
phononic thermal conductivity both in Sr 2 Ru04 
and in «-(ET) 2 salts. The c-axis thermal conduc- 
tivity in Sr 2 Ru04 |H3] appears to be described in 
terms of the phononic thermal conductivity |41j . 
As to k-(ET) 2 Cu(NCS) 2 it appears the phononic 
thermal conductivity dominates for T > 0.5K 

m 



I41| . For T < 0.47K there appears a clear 
sign of the electronic contribution. From the an- 
gle dependence of the thermal conductivity, we 
can deduce d+s-wave (A(k) cx cos(2</>) — 0.067) 
for k-(ET) 2 Cu(NCS) 2 . In Fig.l we show |A(k)| 
for superconductors a)s-wave, (b) d-wave as in 
HTSC, CeCoInggg, k-(ET) 2 saltsjEj, (c)2D f- 
wave as in Sr 2 Ru04 |26J27j and (d)/-wave as in 
UPt 3 |4dl44l45| . 
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2. Thermal conductivity in the vortex 
state. 

In the following we limit ourselves to nodal 
superconductors; 2D /-wave with A(k) oc 
cos(cfc z )e ±I< ^ as in Sr2Ru04 and d + s-wave with 
A(k) oc cos(20) - r as in k-(ET) 2 Cu(NCS) 2 . 
Also we limit ourselves to the clean limit(i.e. 
VAT > e = vVeH, where v = , |57| . Then 
the thermal conductivity for T -C A and in a 
planar magnetic field is given by 
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for <i + s-wave. 

Similarly the Hall thermal conductivity is given 

by 
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for d + s-wave. Here kq is the thermal conductiv- 
ity in the absence of the magnetic field. 

3. Phononic thermal conductivity 

So far we have considered only the electronic 
thermal conductivity. In general, the thermal 
conductivity is written as K = n el + n 9 , where the 
second term is due to phonons. The importance 
of k 9 in high T c cuprate superconductors, YBCO 
and Bi2212 have been well-documented I 11 
low temperatures (for T < 5K ) k 9 ~ T 3 in single 
crystals of YBCO and Bi2212 gHj. In high qual- 
ity crystals phonons are mostly scattered by crys- 
talline defects and crystal boundaries. In the vor- 
tex state in nodal superconductors the quasiparti- 
cles provide another scattering center. When the 
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Figure 2. The field dependence of the thermal 
resistance in the magnetic field along [110] direc- 
tion within the plane Ref. |47| is compared with 
the fitted aH * log(l/b * H) and >/H. 



thermal phonons are ballistic, the phonon scat- 
tering due to the quasiparticles is proportional to 
[N(O,0)] 2 at low temperatures |H]. Here N(0,9) 
is the residual density of states in the presence 
of a magnetic field and 8 refers to the field ori- 
entation. Therefore the c axis phononic thermal 
conductivity k 9 in Sr2Ru04 is a planar magnetic 
field is written 
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for 2D /-wave, and 
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for d + s-wave, respectively. Here To is a con- 
stant of the dimension of the energy, and c is 
the phonon velocity. It is noteworthy that Eq.(8) 
does not contain the fourfold term. 

When the phononic thermal conductivity dom- 
inates as in Si'2Ru04 |47| . it is more convenient 
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T=0.30 K (Suzuki et al.) 
Fitted with a*(H-H cl ) 
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Figure 3. The interplane field dependence of the 
thermal resistance Ref . |47] is compared with a * 
(H-H cl ). 



to analyze the thermal resistance W = (k 9 ) 1 
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for 2D /-wave. We believe aH ln(l /bH) depen- 
dence is more consistent with the experiment for 
H || [100], which indicates the crystal of Sr2Ru04 
is in the clean limit. This is readily seen from 
Fig. 2. Also in the clean limit the specific heat 
behaves similarly A7\. On the other hand for 
H || [001], it appears that the superclean limit 
appears, ie. linear in H as shown in Fig. 3. 

More recently the 6-axis thermal conductivity 
of k-(ET)2Cu(NCS)2 in a planar magnetic field 
is reported (THj. First of all for T > 1 K the 
thermal conductivity decreases with increasing H 
indicating the phonon dominance, more surpris- 
ing is an appreciable twofold term(~ cos(26>)) in 
the thermal conductivity. This suggests naturally 
d + s-wave model. Indeed we obtain A(k) ~ 
cos(20) — 0.067 as already mentioned earlier. As 
the temperature is lowered through T = 0.47 K, 
there appears a positive fourfold term. Clearly 



this is the signature of the electronic thermal con- 
ductivity. Further the magnitude of the fourfold 
term is consistent with Eq.(4). Therefore we con- 
clude the order parameter in k-(ET)2Cu(NCS)2 
is of d + s-wave with A(k) ~ cos(20) — 0.067. In 
other words the nodal lines run 43.08° from the fa- 
axis rather than 45°. Since k-(ET)2 salts do not 
have the tetragonal symmetry, the admixture of 
s-wave component is allowed. This is somewhat 
similar to YBCO, where the tetragonal symmetry 
is broken due to the orthorhombic distortion^]. 
This d x 2_ y 2-w&ve is totally unexpected theoret- 
ically |16I17I18I19| . There is a well-known paral- 
lel between high T c cuprates and /c-(ET)a salts. 
However, the present result tells there is a subtle 
and delicate difference between these two systems 
and a more careful study of the pairing interac- 
tion is necessary. This makes the physics in or- 
ganic superconductors all the more interesting. 

4. Concluding Remarks 

The gap symmetry is the central issue of new 
superconductors. We have shown the angular de- 
pendent thermal conductivity in the vortex state 
provides a new window to look this question. In 
this way we have succeeded in identifying the gap 
symmetry of Sr 2 RuC>4, CeCoIn 5 |42I49| and K- 
(ET) 2 Cu(NCS) 2 H3. We expect this method will 
be very useful to identify the gap symmetry of 
/3-(ET) 2 salts, A-(ET) 2 salts and other organic 
superconductors and to clarify existing contro- 
versy. Also the success of this method testifies 
the soundness of both the BCS theory of nodal 
superconductors and the Volovik's semiclassical 
approach to handle the vortex state in nodal su- 
perconductors. 
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